Role of the trigonal warping on the minimal conductivity of bilayer graphene 
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Using a reformulated Kubo formula we calculate the zero-energy minimal conductivity of bilayer 
graphene taking into account the small but finite trigonal warping. We find that the conductivity 
is independent of the strength of the trigonal warping and it is three times as large as that without 
trigonal warping, and six times larger than that in single layer graphene. Although the trigonal 
warping of the dispersion relation around the valleys in the Brillouin zone is effective only for low 
energy excitations, our result shows that its role cannot be neglected in the zero-energy minimal 
conductivity. 
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Recent experiments have proved that the charge car- 
riers in graphene (single or stacks of atomic layer of 
graphite) are massless Dirac fermion^ii^i^. For recent 
reviews on graphene see Refs. Besides the un- 

usual transport properties observed and reviewed in 
the above works another important experimental fea- 
ture is the minimal conductivity of the graphene sys- 
tems which was considered theoretically^ long before the 
experimental evidence. After the above mentioned ex- 
perimental works on graphene, number of theoretical 
Studies^iaiiaiiia2i3ii4|i5a6^i7ii8 j^^ve predicted the con- 
ductivity of the order of e'^/h. Very recently, Miao 
et al. have experimentally confirmedi^- most theoretical 
prediction !] "'^^d'^d^d^d^d'^d^ , namely the minimum con- 
ductivity in wide and short strips approaches the uni- 
versal value cr^'" = (4/7r) e^/h in single layer graphene. 

The bilayer graphene has been studied first 
experimentally'^ by Novoselov et al. and theoretically^^ 
by McCann and Fal'ko. McCann have calculated the 
asymmetry gap in the electronic band structure of 
bilayer graphene^ ^. In biased bilayer graphene it was 
demonstrated that the gap can be tuned by electric field 
effec1j2^. In bilayer graphene the semiconductor gap has 
recently been controlled experimentally by Ohta et al;^. 
The optical and magneto-optical far infrared properties 
of bilayer graphene has been studied by Abergel et ali^. 
The role of the impurities in biased bilayer graphene 
has been studied by Nilsson and Neto2^. Ludwig has 
considered the conductance of a normal-superconductor 
junction in bilayer graphene'^^-. Recently, Koshino and 
Ando have investigated the transport in bilayer graphene 
in self-consistent Born approximation^! and they found 
that in the strong-disorder regime cr™" = (8/7r)e^//i, 
while in the weak-disorder regime crjf!j,™ — (24/7r)e^/ft- 
which is six times larger than in single layer graphene. 
Similarly, Katsnelson has also calculated the minimal 
conductivity in bilayers using the Landauer approach^^ 
and he obtained a different value a^^" = 2e^/h. In 



Ref . [ITI we found cr'^^^ — (S/tt) e^/h which was confirmed 
later by Snyman and Beenakker— using the Landauer 
approach. 

However, much fewer theoretical works paid attention 
to the role of the trigonal warping in bilayer graphene. 
The influence of the trigonal warping on the weak lo- 
calization effect has been investigated by Kechedzhi et 
alj2£, while on the minimal conductivity only by Koshino 
and Andc^^ using an effective 2 by 2 Hamiltonian. Our 
aim in this work is to calculate the minimal conductivity 
using the Hamiltonian suggested originally by McCann 
and Fal'koSS.. This Hamiltonian allows us to find the 
zero-energy minimal conductivity as a function of the 
strength of the trigonal warping in bilayer graphene. We 
use the Kubo formula rewritten in a form suitable for ob- 
taining the zero-energy minimal conductivity in graphene 
systems. Surprisingly, we find that the conductivity is in- 
dependent of the strength of the trigonal warping and six 
times as large as that for single layer graphene. 

The bilayer graphene consists of two coupled honey- 
comb lattices with basis atoms Ai , Bi and A2 , B2 in the 
bottom and the top layers, respectively. The two lay- 
ers are arranged in Bernal stacking {A2 — Bi). The in- 
tralayer couphng between Ai and Bi, and A2 and B2 is 
7o. The strongest interlayer coupling is between A2 and 
Bi with coupling constant 71 . A direct hopping between 
Al and B2 is taken into account by the coupling constant 
73 ^ 7i- This coupling is responsible for the trigonal 
warping. The above coupling constants are estimated as 
70 = 3.16eVai, 71 = 0.39eV32, and 73 = 0.315eV^. 

To model the bilayer graphene we use the same gap- 
less Hamiltonian as that in Ref. 20, which takes into 
account the trigonal warping. The Hamiltonian in the 
basis Al, Bi, A2, B2 in the valley K and in the basis 
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Bi, Ai, A2 in the valley K' reads 
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where p± = px i ipy, v = \/3a7o/(2fi,) and W3 = 
V^ajz / {2h) , while ^ = +1 for the valley K and ^ = — 1 
for the valley K' (a — 0.246 nm is the lattice constant 
in the honeycomb lattice). The strength of the trigonal 
warping is desribed by the parameter (3 = v^/v = 73/70- 
According to previous studie a^^i'^^'^'^ /3 « 0.1. 

The four eigenvalues of the Hamiltonian ([1]) as func- 
tions of the wave number k = k(cos , soup) are given 

by 

El{k,p) = ^ [1 + P + 2) + (-l)"r] , where (2) 



r 



l-2fc2 (/32 _ 2) + fc4/32 (/32 _,_4)_,_8P^cos3^, 



where n ~ 1,2, while the rescaled wave number is 
k = k^i/{hv). Owing to the cos3iy9 term the eigenvalues 
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FIG. 1: Constant energy lines (in units of 71) of the dispersion 
relation of the positive eigenvalue Ei in the (kx,ky) plane 
around the K point of the Brillouin zone (at the origin in 
this figure). Here /3 = 0.1 and the contour lines are plotted 
equidistantly with the most outer contour line corresponding 
to energy 2El- 

are three-fold rotational invariant for finite [3. The eigen- 
values ±Ei become zero at the K point of the Brillouin 
zone, ie, at fc = 0, and at the center of the three pockets 
located a,i k — (i and p = 0, 27r/3, 47r/3. Around these 
zeros the constant energy lines are distorted as shown 
in Fig. [TJ This is called trigonal warping. At moder- 
ate energy, direct hopping between Ax and leads to 
trigonal warping of the constant energy lines about each 
valley, but at an energy E less than the Lifshitz energy 



El = 71/3^/(4-^/32) « ImeV, the effect of trigonal warp- 
ing is dramatic. It leads to a Lifshitz transitionM: the 
constant energy line is broken into four pockets, which 
can be referred to as one central and three leg parts. 
For V3 ^ w, ie, /3 ^ 1 we find that the separation of 
the 2D Fermi line into four pockets would take place for 
very small carrier densities n < ~ 1 x lO-'^^cm"^. For 
n < riL, the central part of the Fermi surface at energy 
E is approximately circular with area Ac ~ ttE'^/ (Hvs)'^, 
and each leg part is elliptical with area Ae ~ g^c- For 
E <Si El the dispersion relation are linear in fc. A similar 
structure of the constant energy lines can be seen around 
the K' point. For P — there is no trigonal warping, ie, 
the eigenvalues are rotational symmetric and the Dirac 
cones are only at the K and K' points. 

Recently, in self-consistent Born approximation 
Koshino and Ando^"^ have investigated the minimal con- 
ductivity for bilayer graphene using an approximated 2 
by 2 Hamiltonian which mimics the trigonal warping and 
is given by 



Hb2 - 92 ( p2 
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where the effective coupling constant 52 = 7173/70 j the 
rescaled momentums are p± = {px ± ipy)/po and po — 

2/i7i73/(V3a7o)- 

The simplest effective Hamiltonian valid for _E <C 71 
and first introduced by McCann and Fal'kc^l to study 
the Hall conductivity of bilayer graphene is given by 
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where 93 = v'^/ji is the effective coupling constant. In 
this case the trigonal warping is absent. In this work, 
we calculate the minimal conductivity for all the three 
Hamiltonians, Hti,Hi,2 and Hf,^. 

To find the minimal conductivity for graphene systems 
we start from the Kubo formula used by Ryu et al. in 
Ref.[l^ which at zero temperature and for dc conductivity 
(at zero frequency lo) is given by 

o'u!^" = "^sTiv lim cTui/C??), where (5a) 

<y,.Ari) ^ -S,.u^J dVE^,(r,r';i5 = 0,r7),(5b) 

S^,(r,r'; i?,r;) = Tr [G^-i^(r, r'; i?, 77)^^ 

xG^-''{vy:E,r,)j,]. (5c) 

Here (/i, ly) — x, y, the spin degeneracy is Us = 2, the 
valley degeneracy corresponding to the valley K and K' 
is Uy — 2, the area of the sample is S, while 

G^-i^(r, r'; E, 77)^- (r, r'; E, 77)^6+ (r, r'; E, 77). (5d) 

The trace is taken over the spinor indices and for systems 
with translation invariance, the single-particle Green's 
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functions are given by 
G±(ri,r2;i?,77) 
G±(k;i?,ry) 



(2^) 

H{k) = H{p^hk), 
and the current operator is 

.e , e dH{k) 

h 



.k(r.-r,)G±(k;iJ,^),(5e) 
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(5g) 

(5h) 



The above expression (j5bp for the conductivity can be 
simphfied using the identity 

(-Z - ny^ - (z - i?)"^ -2z(z2 - ij2)-\ (6) 

Then with z — iij we can rewrite a^i,{ri) in Eq. (j5b[) as 
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7] I{ri), where 
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This expression of the conductivity is vahd for transla- 
tional invariant systems. 

Before we turn to the case of the bilayer it is instruc- 
tive to see how the expression (j7bp works for single layer 
graphenc. In this case the Hamiltonian is given by 
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where gs = hv and k± = fcj; ± iky. The integrand in 
Eq. (j7b[) can easily be calculated using the polar coordi- 
nates k = fc(cos iy9, sin (p) and one finds 



dk 
^2^ 
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27rr/2 ' 
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which is independent of the coupling constant gs- Note 
that the main contribution in the integral liji) comes 
from the vicinity of fc = 0, therefore the integral over 



k can be extended to infinityi^. Then from Eqs. ([7a|) 
and (jSap we obtain the well-known universal value of the 
minimal conductivity for single layer graphene: crj??^™ = 
(4/^) {eyh)i^^^^^^^^^^2^. 

We now consider the bilayer graphene taking into ac- 
count the effect of the trigonal warping. For bilayer with 
Hamiltonian ([1]) the current operators defined by Eq. (jShp 
have a simple form 



Jx 



/ 1 /3 ■ 

10 

1 

V/3 1 



/ -10/3 

10 

-1 

V-/3 1 



(10a) 



(10b) 



The integral I{rj) = {f3, rj) in Eq. l|7bp will depend on 
(3. It can be shown that a^x — '^yy^ therefore it is conve- 
nient to calculate axx = is'xx + o'ya)/2- Using the polar 
coordinates k = /(;(cos(/?, sintp) and rescaling the vari- 
ables k and rj as k khv/ji and rj — > 77/71, a straight- 
forward algebra yields for the case of ^ = +1 (valley K): 



A + B cos 3ip dip dk 



, (11a) 



/o Jt) {C + DcosSiff 27r 2n' 
A^k'^{2 + 5/3^) 

+ + ry^) [4fc2 (1 + /32) + 2rj^ + P^l + 7?')](Jlb) 
B = 4fc3/3, (11c) 
C = k^ + 7j^ + r]^ + k^ [2jf + Z?' (1 + if)] , (lid) 
D = -2fc3/3. (lie) 

This expression has a three-fold rotational symmetry as 
should be for trigonally warped bilayer graphene. It can 
be shown that for C = ~1 (valley K') we have the same 
results. The conductivity is two-fold degenerate accord- 
ing to the valleys, ie, = 2. 

The integral over ip can be performed analytically and 
Ibi (/3, 77) becomes 

./ON ^ r , AC - BD dk 

IbiiP,T]) = 2/ k — — . 12 

Jo (C2 - D2)3/2 27r 

Without trigonal warping, ie, for (3 — one finds 



4i(/3 = 0,77) = 
dk ^k (fc-* -t- T] 
2^ 



2fc2 + 2fc2,y2 



{k'^ + r]'^ + rf- + 2k'^r]'^y 



(13) 



Thus, using Eqs. (I7ap and (|5ap the minimal conduc- 
tivity for bilayer graphene without trigonal warping is 
crJ"^"'(/3 = 0) = (8/^7r) (eV/i). This result has been de- 
rived first in Rcf. 17 in a different way, and subsequently 
by Snyman and Bcenakker in Ref. 2^ using the Landauer 
approach. 

After a tedious calculation the integral in Eq. for 
finite value of /? can be performed yielding 
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Again, using Eqs. (j7ap and (|5ap we find a remark- 
able result, namely the minimal conductivity for bilayer 
graphene with trigonal warping takes a universal value 
cr^^"'(/3) = (24/7r) {e^/h) independent of the stregth /3 of 
the warping. This is our central result in this paper. This 
value is six times as large as the conductivity in single 
layer graphene. It is suprising that (t™"(/3) is not a con- 
tinous function around /3 = 0. Indeed, as we have seen 
cr™"(/3 = 0) = (S/tt) (e2//i), while for any finite values 
of (3 it is three times larger. This non-analitic behaviour 
of <J^™i(3) at /3 = is a consequence of the fact that the 
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minimal conductivity results from the electronic dynam- 
ics in the limit of zero density n —^ 0. For any non-zero /3, 
such density is always below the Lifshitz density n < ul 
where the 2D Fermi line around each valley forms four 
separate pockets, whereas for /3 — 0, the Lifshitz tran- 
sition does not occur and there is always a single Fermi 
line at each valley. 

In the frame work of self-consistent Born approxi- 
mation the same result was predicted by Koshino and 
Andc2I using the Hamiltonian iJ{,2 given by Eq. ^ . Note 
that in this Hamiltonian there is no adjustable parameter 
for the strength of the trigonal warping like (3 for Hamil- 
tonian ID). The effective coupling constant g2 drops out 
in Eq. (|7bp . therefore in this model the trigonal warping 
is built in without the possibility to change its strength. 
Using Eq. ([7]) we repeat the calculation with the Hamil- 
tonian Q and find 
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(12 - 1277]^) + 0(lnr/) + 0(t]^). 

(15) 



Thus, the minimal conductivity takes the same universal 
value cr5f^"(/3) = (24/7r) (e^/h) as that for the Hamilto- 
nian Hbi given by 

Finally, we calculate the minimal conductivity using 



the simplest Hamiltonian iJ^a given by Eq. (|4]). Then 
the integral in (|7b[) becomes 



2l3 



dk 8.93 fc 



27r 772 -I- gjk-^ 



= A- (16) 



Thus, the minimal conductivity takes the same universal 
value CT™" = (S/tt) (e^/ft.) as that for Hamiltonian ([T|) 
with P = 0. 

In summary, we compared the minimal conductivity in 
bilayer graphene obtained from three different effective 
Hamiltonians used in the literature. We found that for 
the case when the trigonal warping is absent, the conduc- 
tivity is always two times larger, while in the presence of 
trigonal warping it is six times larger than that for single 
layer graphene and is independent of the strength of the 
warping . Our universal results suggests that the con- 
ductivity has a topological origin, which can be a further 
research topic in the future. 

We gratefully acknowledge discussions with C. W. J. 
Beenakker, E. McCann, and V. Fal'ko. This work is 
supported partially by European Commission Contract 
No. MRTN-CT-2003-504574 and the Hungarian Science 
Foundation OTKA 046129. 



^ K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. 

Zhang, S. V. Dubonos, I. V. Grigorieva, and A. A. Firsov, 

Science 306, 666 (2004); K. S. Novoselov, A. K. Geim, S. 

V. Morozov, D. Jiang, M. I. Katsnelson, I. V. Grigorieva, 

S. V. Dubonos, and A. A. Firsov, Nature 438, 197 (2005). 
^ Y. Zhang, J. P. Small, M. E. S. Amori, and P. Kim, Phys. 

Rev. Lett. 94, 176803 (2005); Y. Zhang, Y.-W. Tan, H. L. 

Stormer, and P. Kim, Nature 438, 201 (2005). 
3 K. S. Novoselov, E. McCann, S. V. Morozov, V. I. Fal'ko, 

M. I. Katsnelson, U. Zeitler, D. Jiang, F. Schedin, A. K. 

Geim, Nature Physics 2, 177 (2006). 

M. I. Katsnelson, Materials Today 10 20 (2007). 

^ M. L Katsnelson and K. S. Novoselov, |cond-mat/0703374| 
® A. K. Geim and K. S. Novoselov, Nature Materials 6, 183 

(2007). 

E. Fradkin, Phys. Rev. Rev. B 63, 3263 (1986); P. A. Lee, 
Phys. Rev. Lett. 71, 1887 (1993); E. V. Gorbar, V. P. 
Gusynin, V. A. Miransky, and I. A. Shovkovy, Phys. Rev. 
B 66, 045108 (2002). 

* K. Ziegler, Phys. Rev. Lett. 97, 2668 02 (2006). 

® L. A. Falkovsky and A. A. Varlamov, "cond-mat /0606800' 
^° K. Nomura and A. H. MacDonald, cond-mat/0606589 
" V. P. Gusynin and S. G. Sharapov, Phys. Rev. LettT 95, 
146801 (2005). 

^2 V. P. Gusynin and S. G. Sharapov, Phys. Rev. B 73, 
245411 (2006). 

N. M. R. Peres, F. Guinea, and A. H. Castro Neto, Phys. 
Rev. B 73, 125411 (2006). 

M. I. Katsnelson, Eur. Phys. J. B 51, 157 (2006). 

P. M. Ostrovsky, I. V. Gornyi, and A. D. Mirlin, Phys. 

Rev. B 74, 235443 (2006). 

J. Tworzydlo, B. Trauzettel, M. Titov, A. Rycerz, C.W.J. 



Beenakker, Phys. Rev. Lett. 96, 246802 (2006). 

J. Cserti, Phys. Rev. B 75, 033405 (2007). 

S. Ryu, C. Mudry, A. Furusaki, and A. W. W. Ludwig, 

cond-mat/0610598 

F. Miao, S. W ijeratne, U. Coskun, Y. Zhang, and C. N. 
Lau, |cond-m at /0703052 

E. McCann and V. I. Fal'ko, Phys. Rev. Lett. 96, 086805 
(2006). 

E. McCann, Phys. Rev. B 74, 161403 (2006). 
E, V. Castro, et al. cond-mat/0611342 
T. Ohta, A. Bostwick, T. Seyller, K.^Horn, and E. Roten- 
berg. Science 313, 951 (2006). 



D. S. L. Abergel and V. I. Fal'ko, [cond-mat/0610673 
J. Nilsson and A. H. Castro Neto, |cond-mat/0612'470 
T. Ludwig, cond-mat/0701045 



M. Koshino and T. Ando, Phys. Rev. B 73, 245403 (2006). 
M. L Katsnelson, Eur. Phys. J. B 52, 151 (2006). 
I. Snyman and C. W. J. Beenakker, Phys. Rev. B 75, 
045322 (2007). 

^° K. Kechedzhi, V. I. Fal'k o, E. McCann, and B. L. Alt- 
shuler, cond-mat/0701690 

W. W. Toy, M. S. Dresselhaus and G. Dresselhaus, Phys. 
Rev. 109, 272 (1958). 

A. Misu, E. Mendez, and M. S. Dresselhaus, J. Phys. Soc. 
Jpn. 47, 199 (1979). 

R. E. Doezema, W. R. Datars, H. Schaber, and A. Van 
Schyndel, Phys. Rev. B 19, 4224 (1979). 
A. A. Abrikosov, Fundamentals of the Theory of Metals, 
Elsevier Science Publishers B. V., North-Holland, Amster- 
dam, 1988. 



